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Abstract 

The endochronic theory, developed in the early seventies, allows the plastic behavior 
of materials to be represented by introducing the notion of intrinsic time. With dif- 
ferent viewpoints, several authors discussed the relationship between this theory and 
the classical theory of plasticity. Two major differences are the presence of plastic 
strains during unloading phases and the absence of an elastic domain. Later, the en- 
dochronic plasticity theory was modified in order to introduce the effect of damage. 
In the present paper, a basic endochronic model with isotropic damage is formu- 
lated starting from the postulate of strain equivalence. Unlike the previous similar 
analyses, in this presentation the formal tools chosen to formulate the model are 
those of convex analysis, often used in classical plasticity: namely pseudo-potentials, 
indicator functions, sub-differentials, etc. As a result, the notion of loading surface 
for an endochronic model of plasticity with damage is investigated and an insight- 
ful comparison with classical models is made possible. A damage pseudo-potential 
definition allowing a very general damage evolution is given. 
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Introduction 



In the early seventies, IValanisI (119711 ) proposed the endochronic theory of visco-plasticity, 
which postulates the existence of an intrinsic time governing the rate-independent evolu- 
tion of stress and strains in materials, whereas the Newtonian t ime i s exploited to model 



the viscous behavior; see also (ISchapery 



1968 



Bazant and Bath 



19761 ) . In the case of plas- 



ticity without viscous effects, the resulting constitutive laws are characterized by the ab- 
sence of an elastic domain and the corresponding hysteresis loops are typically smooth and 
open. The flow rules of these models were not originally formulated in terms of pseudo- 
potentials, which ma de the direct comparison of this class of models wi th classical plastic- 



i ty th eories difficult (IValanis 



Erlicher and Point 



19801 ) . However, it was recently proven by 
(J2006|) that endochronic models do admit a representation based on pseudo-potentials and 
on the normality assumption, provided that pseudo-potentials be endowed with an additional 
dependence on state variables. This proof, given for the case of plastically incompressible 
models, showe d the strong relationship between the endochronic th e ory and the general- 



ized plasticity (jPhillips and Sierakowskil . 



1993 



Auricchio and Taylor! . 



1965 



Eisenberg and Phillipsl . 



1971 



Lubliner et al. 



19951 ). It was also shown that the non-linear kinematic hardening 



model, that is associated, but is not in a generalized sense, admits a representation in terms 



of a pseudo-potentia 
like the Mroz mode 



(lErlicher and Point 



. Recently, the same au thors extended this analysis to other models, 



Point and Erl 



20051 1: see also 



cherl. 



2007) and the non- a s sociated Drucker-Prager m odel 



Ziegler and Wehrlil fll987f ). 



Houlsby and Puzrinl (120001 ) . In 



summary, this thermodynamically well-posed approach can be used for a very large class of 
existing classical or non-classical plasticity models. Actually a similar approach is used in 



geotechnical engineering, see e.g. 



Collins and Houlsbyl (119971 ). where pseudo-potentials have 



an additional dependence on the so-called true stresses, distinguished from the generalized 
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stresses. 



The standard endochronic theory was modified by several auth ors thro u gh th e introduction 



of a damage variable. Using the strain equivalen ce postu l ate, 



endochronic model with isotropic damage, while IValanid (119901 ) discussed an endochronic 



Xiao del (119891 ) proposed an 



model with anisotropic damage, in the larger theoretical framework of fracture mechanics. 



Later, a d ifferent approach based on the postulate of en e rgy equiva 



others, by I Chow and Chenl (119921 ) and 



Wu and Nanakornl (11998 



19991). 



ence was used, among 



In the aforementioned works, the thermodynamic formulation of flow rules is not based on 
the notions of pseudo-potentials and loading surfaces, as it is typical for other classical plas- 



ticity models with or without damage. Hence, in this paper, a simp 



e endochronic model of 



Xiao del (119891 ) is presented: no 



plasticity with isotropic damage similar to that discussed by 
generalization is introduced with respect to the previously cited models, but a new approach 
is suggested for their description. In detail, the postulate of strain equivalence is adopted; the 
Helmholtz energy is assumed to have a regular quadratic term and an additional si ngular term ; 



the tools of th 


ie convex analysis such 


1969; 


Moreau 


1970; 


Fremond . 


2002) 



20021 ) are used to define the flow rules starting from well-suited 
pseudo-potentials. This presentation leads to the proper definition of the plasticity loading 
surface for an endochronic model with damage and is a direct exte nsion of the results concern - 



Erlicher and Point 



(|2006j). 



ing the endochronic model without damage already discussed in 
Only plastically incompressible models are considered here, since they permit to explain the 
main ideas, without introducing a too complex formalism. The extension to the general case is 
possible, but it is omitted for simplicity. The proposed analysis has an intrinsic interest, since 
it allows an easier comparison between endochronic models with damage and classical plas- 
ticity models with damage. Nonetheless, in the authors' opinion, another important reason 
justifies the interest towards this class of models: they represent the suitable theoretical basis 
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for the analysis of the thermodynamic admi ssibility of the Bouc-Wen models with st r ength 



and stiffness degradation; see among ot hers (IBou 



Casciatil . 



1989 



1971 



Wen 



1976 



Baber and Wen 



1981 



Karrav and Bouc 



19891 ) . This was one of the main motivation at the origin 



of the present study and the rel ated developments about degrading Bouc-Wen models are 



presented in a companion paper (jErlicher and Bursi 



2007J) 



After the introduction, the endochronic theory is presented in the second section: in the first 
part, standard endochronic models are described, while the second part concerns the definition 
of the flow rules of the extended endochronic theory, characterized by an additional scalar vari- 
able endowed with damage. The thermodynamic framework, with the definition of the suited 
pseudo-potentials, is discussed in the following section and is supplemented by numerical ex- 
amples. Then, a brief discussion about stability and uniqueness is made and the concluding 



rema rks are given, where the topics dealt with in the companion paper (jErlicher and Bursi 



20071 ) are pointed out. 



Endochronic models 



Flow rules of plastically incompressible ND-EC models 



Valanid (I197ll ). who suggested the use of a 



The endochronic theory was first formulated by 
positive scalar variable called the intrinsic time scale, in the definition of constitutive plas- 
ticity models. The evolution laws are described by convolution integrals involving past values 
of the strain e and suitable scalar functions /i depending on called memory kernels. When 
the memory kernel is exponential, the integral expressions can be rewritten as simple differen- 
tial equations, the flow rules; in the case of an isotropic endochronic model without hardening 
or softening, called here ND-EC model (see Figure 1), fulfilling the plastic incompressibility 
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assumption, they read: 



tr (cr) = 3K tr (e) , dev (cr) = z 
z = 2G dev (i) -(3 z$ 



(1) 



where (3 > (notice that /3 different from zero is needed to have a non elastic behavior); the 
superposed dot indicates the time derivative; e is the small strain tensor; cr is the Cauchy 
stress tensor; tr and dev are the trace and deviatoric operators; K is the bulk modulus while 
G is the shear modulus. The simplest choice for the intrinsic time scale flow indicated in (1) 
is •& = \\dev It is interesting to note that relationships (1) are equivalent to 



cr = C : (e - e p ) 



(2) 



tr (i p ) = 



and 



2G 



where the trace of the plastic strain flow e p is zero, consistently with the assumption of 
plastic incompressibility. C = (K — 2G/3) 1 <g> 1 + 2GI is the elasticity fourth-order tensor for 
isotropic materials; 1 is the second-order identity tensor; I is the fourth-order identity tensor 
and <S> represents the tensor product. 

Flow rules of plastically incompressible D-EC and DD-EC models 



An endochronic model with isotropic hardening or softening with plastically incompressible 
flow is defined as follows: 



tr (cr ) = 3K tr (e ) , dev (cr ) = z 

z = 2G dev (e) - (3 z ■& with & = j 



(3) 
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where g > is called the hardening- softening junction (jBazantl . Il978l ). As stated by its name, 
the function g introduces isotropic hardening (or softening), which distinguishes this model 
(D-EC) from the basic ones presented in the previous section and indicated as ND-EC (see 
Figure 1). In the classical endochronic formulations, g is a functi on of C, where ( is the intrinsic 



Valanid ( 1l97ll ). Another more 



time measure. A standard choice is ( = \\dev according to 

general definition, leading to a cyclic behavior similar to that of the Prandtl-Reuss model 



(jLemaitre and Chaboche 



19901 ) when the positive parameter n is large enough, reads 



c 



1 + — sgn (z : dev {i))J |z : dev (e)| ||z|| n 2 



(4) 



with 7 G [— (3, f3\ in order to ensure the non-negativity of (; sgn is the signum function. An 
important difference between (4) and the standard definition ( = \\dev(i)\\ is related to the 
product z : dev (i), entailing ( = when the deviatoric strain increment is orthogonal to the 
stress. However, ( can be different from zero during unloading, i.e. when z : dev (e) < 0. Eq. 
(4) shows that 7 affects the difference between the loading and unloading values of the intrinsic 
time increment at a given stress z. In particular, when 7 = /3 these increments are zero during 
unloading, while 7 close to (and greater than) —f3 leads to relatively small increments during 
loading, while ( is relatively large during unloading. The influence of n on the endochronic 
model behavior is discussed in the last Section, with reference to the strain accumulation and 
the stress relaxation effects. According to (3) and (4) and assuming (3 + 7 > 0, the norm of 
the tensor z (t) is bounded as follows: 

' 2G \» 



\dev (a (t))\\ < a u 



(5) 



for t > 0, provided that ||z(0)|| < o u . This inequality proves that a limit strength value 
exists and only concerns the deviatoric part of the stress cr, consistently with the plastic 
incompressibility requirement. Eq. (5) also shows that this bounding stress depends on the 
parameters (3, 7 and n. 
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The expression (3) is equivalent to 



a = C : (e - e p ) 

tr (i p ) = and 



(6) 



ip dev(a) £ 

' _ 2G//3 g 



From the last relationship in (6), it appears that the parameters 7 and n introduced in (4) 
affect the amplitude of the plastic strain flow, while the direction is always that of dev (<r). 

A larger class of endochronic models can be defined by the following relationships 



a = (l-D)C :(e-e p ) 
tr(i p )=0 and 



(7) 



1 dev(cr) £ 
l-D 2G//3 g 



where D is a scalar variable introducing isotropic damage. The plasticity model with damage 
defined by (7) is named here the extended endochronic model and it belongs to the class 
of DD-EC models, as indicated in Figure 1. Note that the stress is defined by introducing 
the factor (1 — D), consistently with the d e finitio n of effective stress and the principle of 



strain equivalence ( ILemaitre and Chabochei . 



1990). Moreover, it can be observed that the 



relationships (7) are equivalent to 

tr (tr) = (1 - D) 3K tr (e) , dev (a) = z 
z = (1 -D) 2G dev(i) -(3 z ■& - D ^ with ■& = £ 

which can be compared with (3). 

A possible choice for <^ is given by 



(8) 



C = ( 1 + ~j s 9 n ( z : dev (e)) j |z : dev (e)| ||z|| n 2 (1 — D) 1 



(9) 



S 



which represents a direct generalization of (4): the last factor depending on D and n is 
introduced in order to have an intrinsic time depending on the effective stress instead of the 
actual one, consistently with the strain equivalence postulate. An elastic with damage model 
can be defined by assuming ( = 0. In the authors' knowledge, the notions of pseudo-potential 
and loading surface were never applied to the extended endochronic theory; therefore, these 
aspects are analyzed in detail in the next section. 



A thermodynamic framework for the extended endochronic theory 

The aim of this section is to define the Helmholtz free energy and the pseudo-potential leading 
to the flow rules (7) or, equivalently, (8). Under the assumption of isothermal and small 
transformations, the Helmholtz free energy density is chosen as follows: 



ifr = * ( v ) = 4 ( e> £ P ,(,D) + l n (e, e p ,(,D) 



(10) 



where v = (e, e p ,(,D) is the vector of state variables; e, e p and D were previously defined; 
( is a scalar internal variable associated with isotropic hardening. For all the state variables, 
an initial zero value is assumed. The choice of ( to indicate an internal variable might seem 
misleading, since the symbol ( was also used in (3)- (8) to define the intrinsic time measure. 
However, as it will be seen hereafter, this choice is the proper one, as for endochronic models, 
£ has simultaneously both meanings; ip is the regular part of the Helmholtz energy; Ih is 
the indicator function of the closed set EL by definiti on, an indicator function is equal to 



inside H and equal to +oo outside (IRockafellar 



19691 ); the set EI indicates the admissibility 



domain for the state variables v and should be introduced every time some conditions on 



state variables are to be imposed: for instance, it is equa 



to impose the admissible values for the damage variable (IFremond 



to the interva l D £ [0, 1] in order 



20021 ). 
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Once \I/ is known, the non-dissipative thermodynamic forces q nd = (er nd , T nd , R nd , Y nd ) are 
defined as the gradient of ifj (e, £ P ,(,D): 



nd 



dip 



nd 



S± R nd . = 9± ynd . = ^± 

deP' ' dC' '3D 



(11) 



while the non-dissipative reaction forces q' 



ndr 



r ndr T ndr ^ R ndr ^yndr^ ^ giyen by 



(a ndr , r ndr , R ndr , Y ndr ) e dl m (e, e p ,(,D) 



(12) 



where d is the sub- differential operator (jRockafellarl . 



19691 ). If the constraints imposed by H 



are fulfilled, the indicator function Ijg is zero and ^ (t) — tp (t). This entails the identity of 
the time-derivative s, viz. \l/ ( t ) = iJj (t) + q ndr ■ v = ip (t). In other words, one has q ndr • v =0 



for every instant t ( jFremond 



2002). 



Due to the assumptions of isothermal and small transformations, the expression of the second 
principle reads: 

= cr : i-ip > (13) 

(13) states that the intrinsic (or mechanical) dissipation $x has to be non-negative. Intro- 
ducing the dissipative thermodynamic forces q d = (cr d , r d , _R d , Y d ) as 



er := <x — <r — <r 



T := — r — T 



^j>d R R ndr yd ynd yr 



and substituting (14) in (13), one obtains: 



$i (t) = <T d : i + T d : i p + R d ( + Y d D > 



(14) 



(15) 



In order to fulfill (15), the flows of the state variables i,i p ,( and D have to be suitably 
correlated with the dissipative thermodynamic forces <j d , T d , R d and Y d . Therefore, some 
additional complementarity rules need to be defined: usually, a scalar non-negative function 
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called pseudo-potential 



(v'; v; p) such that 



(0; v; p) = for all v and p 



(16) 



is introduced and the dissipative forces q d = (er d , T d , R d , Y d ) are derived imposing the so- 
called generalized normality assumption on it. Equivalently, one can define the flow rules v 
by imposing the generalized normality assumption on th e dual pseudo-potential <f)*, which is 



the Legendre-Fenchel transform of (IRockafellar 



19691 ). This last method will be explicitly 



exploited herein. The generic flow v' is noted with "prime" , while the actual flow at the present 
state is noted with v. As a matter of fact, the pseudo-potential is assumed to vary with the 
present value of state variables v and with some additional parameters collected in the vector 



p. Th ese parameters may be any quantity related to the past history of the material (jFremondl . 
20021 ) . For instance, one may have p(x) = (e(x), ||e(x)|| ) = (e(x), max < t /< t ||e (x, t')\\), 



where e is the dissipated energy per unit volume at the point x of the body volume and ||£|| max 
is the maximum (from t' = to the present state t' = t) of the strain norm at the same point. 
Observe that the parameters collected in p could also be non-local, like p (x) = Jvb(x) e( ^V , i-e. 
the energy dissipated in a given volume Vo (x) around the point x of body volume. 

When no viscous effect occurs, the case of plasticity with damage is recovered. This corre- 
sponds to choose a pseudo-potential 4> independ ent from s, entailing cr d = 0; for a detailed 
derivation of these relationships, see, for instance, 



Erlicher and Point 



(120061 ). Moreover, "plas- 



tic flow may occur without damage and dam age may occur without appreciable macroscopic 



plastic flow" (ILemaitre and Chaboche 



19901 ) . Therefore, (15) with a d = "must be split in 



two independent inequalities": 



e p := T d :i p + R d (> 0, 



e D := Y d D > 



(17) 



The two scalar quantities e„ and e^ respectively define the rate of energy per unit volume 
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dissipated by plasticity-related phenomena and by damage phenomena; see Figure 2. Their 
sum e = e p + en is the rate of the total dissipated energy per unit volume and coincides 
with the intrinsic dissipation $1. The restrictions imposed by these two inequalities are more 
severe than the original unique inequality of Clausius-Duhem (15). However, they are usually 
adop ted as basic thermodynamic cr iterion for the formulation of plasticity models with dam- 



age (iLemaitre and Chaboche 



19901 ). This assumption will be adopted hereafter. Taking into 
account (17), the pseudo-potential is supposed to split into two pseudo-potentials 4>d and (j) p , 
respectively related to damage and plastic flow: 



(>', C', D'\ v; p) = <j) D {if- v; p) + P (e p \ C'; v; p 



;is) 



In the following sections, the Helmhol tz free energy, the 



as their Legendre-Fenchel transforms (jRockafellar 
model with damage (7). 



pseudo-potentials 4>d and <f> p , as well 



19691 ). are formulated for the endochronic 



The Helmholtz free energy 



According to (10), for the DD-EC models one has the following Helmholtz free energy: 



* (v) = V (v) + I H = (1 - D) - (e - e p ) : C : (s - e p ) + h 



In this paper, two cases are considered: 



H = < 



(e, e p , C, D) such that D > 0, D < 1 
and (1 - D) s R (v, p) - r < 



(19) 



(20) 



where s and ro are positive parameters; v = (e, e p , () collects all state variables except D 
and R = R(~v, p) is a non-negative function called source of damage. The first two conditions 
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on D impose the minimum and the maximum values for this variable. As it will be seen, the 
third condition in (20) is strictly related to the definition of the damage limit surface. The 
second case is characterized by a different assumption: 

H = {(e, e p , C, D) such that D > 0, D < 1} (21) 

where only the two basic inequalities on D are retained. 

Making use of (14), (19), (20) and of the pseudo-potential (25), i.e. the Definition 1 of <pD given 
in the following section, it is possible to prove that the assumption q ndr = is admissible. 
The same holds for the model defined by (19), (21) and (34) (Definition 2 of 4>d)- For brevity, 
the details of this proof, are omitted. As a result, the non-dissipative thermodynamic forces 
fulfill the following relationships: 

a nd = || = (1 - D) C:(e-e p )=(T-(T d = (T 



T nd = = — (1 — D) C:(e-e p ) = -t° 



R nd = M = o = -R d 

ynd = ^t = _l ( £ _ £ P) . C . ( £ _ e p) = _yd 

Moreover using (22) and supposing D < 1, the energy dissipation rate reads 



e = e p + e D = a : i p + \a : f~iD 

= dev (a) : dev {&) + ^tr (i p ) + \ ( ^Mg^ +Mgl!) * 



(22) 



(23) 



Taking the time- derivative of (22)i and recalling that no viscous effect is considered (cr d = 0), 
one obtains 

i = iP + i e = i P + c -i. + c -i . __Z := iP + e e > a + i e ' D (24) 

l-D (l-D) 2 V 1 
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where e e ' a is an elastic strain flow at constant damage and i e ' is an elastic strain flow at 
constant stress. It follows that eo = 1/2 (& : i e ' D ^j. Note that in general the endochronic 
theory may present non-zero energy rates e p and en also during unloading phases; see Figure 
2b in this respect. 



The pseudo-potential for the damage flow 



The formalism of the loading function fp, as well as the pseudo-potentia 



be used to express the 



Nedjar . 



2001 



Fremond 



d amag e evolution (ILemaitre and Chaboche . 



1990 



I 



Can 



Salari et al 



2004 : 



20021 ). We present herein a well-known damage evolution rule by 



using both pseudo-potentials 4>r> and its dual <f>* D . Then, a discussion is done about a novel 
pseudo-potential leading to a damage evolution where D may be different from also during 
unloading phases. In detail, the main difference between the two cases is related to the role of 
the damage limit surface. Standard damage evolution rules, viz. Definition 1, are characterized 
by the possibility for the actual state point to be inside the damage domain delimited by this 
limit surface; in this situation and in particular during unloading phases, damage increments 
are null. Conversely, in the damage evolution which we propose here, i.e. Definition 2, the 
present state point is forced to be always on the damage limit surface also during unloading 
phases. 

Definition 1 of 4>d 

Let us begin with the following pseudo-potential, associated with the Helmholtz free energy 
(19)-(20): 



{e - e p ) : C : (e - e p ) - (1 - D) s R (v, p) + r 



D' + I Bd (D') (25) 
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and 



P D = [£>' such that D' > o} 



(26) 



The pseudo-potential <fir> is the sum of a regular part, proportional to and of the indicator 
function Io D . The term multiplying D' in the regular part of 0^ is always non- negative, by 
virtue of the third condition defining EI in (20). The regular part of 4>d, considered for the 
actual flow D' = D, represents the rate of dissipated energy eo- (25)-(26) allow a large number 
of standard damage evolution rules to be represented, according to the specific definition of 
R (v, p). An interesting example is 



R — R(e) — - [2Ge e > + : e 



: > + + A([tr(s e )] + ) 2 ) 



(27) 



where A = K — 2G/3 is the Lame constant. For a scalar x, [x} + := (x), where () are the 
McCauley brackets. The positive part e e,+ of the tensor e e i s obtain e d afte r dia gonalisation. 



Other definitions for R can be adopted; see, among others, 



fl2004h . 



Nediarl (J200l|) and 



Salari et al. 



In order to derive the damage flow, it is convenient to consider the Legendre-Fenchel transform 
of (po, which reads: 



f D (F d ';v;p) =^Vd' &d [( 



I Ed (Y d ';v;p 



Y d ' - \ (e - e p ) : C : (e - e p ) + (1 - D) s R (v, p) - r ) D' 



(28) 

where Eo= ^Y d ' such that fjj (Y d '; v; pj < o| is the damage loading domain and the cor- 
responding loading function is: 



* 1 



f D (Y d '; v; p) = Y d ' - - (s - e?) : C : (s - ^) + (1 - D) s R (v, p) - r := Y d ' - Y± x (29) 
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By using the normality assumption, the relevant damage flow rule reads 



17 — A D d y d > — — *D 



(30) 



with f D (Y d ; v; p) A D = 0, A D > 0, / D (Y d ; v; p) < 0. 

At the actual state, it holds Y d ' = Y d = -Y nd , and therefore: 

Id = fo (Y d ; v; p) = (1 - D) s i? (v, p) - r < 0, (31) 

which is the damage limit surface, but also is one of the conditions defining the set H. It 
becomes evident that the positive constant tq is the initial damage threshold. The Kuhn- 
Tucker conditions state that f D < implies no damage increment, while fo (Y d ; v; pj = 
corresponds to a damage increment which can be computed by enforcing the consistency 
condition: 

j D = (1 - D f (^£> i+ d -^p) - R (*, P) . (1 - Z))'- 1 i» = (32) 
leading to the explicit expression of the damage flow 

where H is the Heaviside function. The presence of the Heaviside function in the damage flow 
definition indicates that damage increments are zero during unloading phases. Note that (33) 
entails that the limit condition D = 1 is never reached. 

Figure 3 illustrates some loading-unloading cycles of an elastic with damage model (( = 0). 
The uniaxial stress is considered, viz. all the components of the Cauchy tensor are supposed 
to be null, except an. The parameter values represent a hypothetical material for which the 
Young modulus E = 35000 MPa and the Poisson ratio v — 0.18 are close to those of concrete; 
damage is defined by (20) and (27), with s = 2.5 and r = 1.2e — 05 MJ/m 3 ; see also the 
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numerical examples in 



Nedjarl (120011 ) . Together with the stress-strain and damage evolution of 



this model, Figure 3 depicts the evolution of Y d = 1/2 (e — e p ) : C : (e — e p ), i.e. the actual 
value of Y d \ and of the quantity Y*^ = 1/2 (e - £ p ) : C : (e - e p ) + r - (1 - D) S R (e), 
defining the upper limit of Y d ' according to (29). When these two curves are superposed, the 
damage increases. 



Definition 2 of 4>d 

Unfortunately, a definition of D of the type (33), deriving from the pseudo-potential (25) 
and the condition (20), is not able to represent the case of damage increasing during both 
loading and unloading phases, owing to the condition f D < 0. We recall that the case of 
damage increasing during unloading may occur in Bouc-Wen models with stiffness degradation 



(jErlicher and Bursi 



20071 ). A damage pseudo-potential, simpler than (25), is more suited: 



</> D (!>'; v) = [i (e - e p ) : C : (e - e p )] D' + I Bd (£>') (34) 

with 3 D still provided by (26) and with the conditions on the damage state variable defined 
in (21). As already observed, it is possible to prove that the assumption q ndr = is admissible 
also for this Definition 2 of the damage pseudo-potential. The dual pseudo-potential becomes 
(jf D (Y d '; v) = I Ed (Y d '] vj where E D = |V d ' such that f D (Y d '; < j is the correspond- 
ing damage loading domain, with the damage loading function 

f D (Y d '; v) = Y d ' - \ (s - e p ) : C : (s - e p ) := Y d ' - Y± x (35) 

At the actual state, Y d ' = Y d = —Y nd and therefore fu (Y d ; vj = at every instant. 
Therefore, the relationships (30) reduce to D — \ D dfi)/dY d ' — \ D , with X D > 0. Moreover, 
D = Xo can no longer be computed by the consistency condition, fulfilled as an identity at 
every instant. Hence, it must be rather defined by an additional condition. Any definition 
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ensuring rate-independence, consistent with (21) and fulfilling D = Ad > is admissible, 
even though is characterized by non-zero damage increments during unloading phases. 



The pseudo-potential for the plastic flow 



The usual method to define associated plastic flows is based on the notion of loading function, 
indicated here by f p , as well as on the normality assumpti on. Another eq uivalent formalism is 



based on the use of the dual pseudo-potential (j>* = I/ p <o ( jMoreaul . 



19701 ) . A third way to for- 



mu late plasticity models is based on the pseudo-potential cj) p , Legendre-Fenchel conjugate o: 



remond 



2002 



Ziegler and Wehrfil 



1987 



Houlsby and Puzrin 



2000 



Erlicher and Point 



20061 ). The advantage of using the formalism based on f p (or 0* = 1/ <o) is essentially simplic- 



ity. Moreover, when a non-associated flow is to be defined, the simple introduction of a second 
function g p called plastic potential matches this purpose. Nonetheles s, for some non-class ical 



plasticity theories, like endochronic theory and generalized plasticity (jLubliner et al 



199.Tl.it 



is not straightforwarc 



by 



Erlicher and Point 



to pro vide a proper definition of the loading function f p . It was proved 



(120061 ) that for these plasticity theories (without damage) a way to de- 
fine the loading function is to start from the definition of the pseudo-potential <j) p , to compute 
the dual potential 0* and then to derive f p . An important point is the additional dependence 
of 4> P , and therefore of 0* and the loading function too, on the state variables. This dependence 
is only optional for standard plasticity theories but is essential both for the endochronic theory 
and the generalized plasticity. Moreover, we notice that some models with non-associated flow 
also admit a representation based on the definition of a suited pseudo-potential <fr p , depending 



on state variables. 



Erlicher and Point 



he e xample of a non-associated Drucker-Prager model can be found in 



(120051 ): in particular, it is shown that a suited pseudo-potential <p p leads 



to a modified loading function which plays both roles of the traditional loading function and 
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of the plastic potential. 



For the endochronic models with damage, the plasticity pseudo-potential is defined as follows: 



(e^ftvjp ) = (!-£>) 



\dev {C : {e - e p ))f (' 



2G/P 



g{y,p 



+ i d (^',C';v;p) (36) 



where In is the indicator function of the convex set 



D = < 



(e p ',C') such that 

tr (>') = 0, (' > and 



, _ dev(C:(e - e p )) 



2G//3 



fl(v,p) 



(37) 



(see Figure 4a). The first equality in D imposes the plastic incompressibility of the flow. 
Moreover, since D is supposed to be less or equal to one and g = g (v, p), the hardening- 
softening function, is positive by assumption, the second condition in D ensures the positivity 
of (j> p . Therefore, the standard properties of <f> p , viz. non-negativity, convexity and positive 
homogeneity of order 1, are matched. The third condition in D gives the plastic flow and is 
consistent with (7). It can be proven that when i p = £ p and (' = C, i-e. when the actual flows 
are considered, the first term of the sum in (36) represents the rate of energy e p dissipated by 
the plastic flow, defined in (17) for the general case. Note that the pseudo-potential has an 
additional dependence on the state variables and on the past-history dependent parameters 
collected in p. 



he dual dissipati on potential </>* is obtained by the Legendre-Fenchel transformation of 



( jRockafellar . 



1969T ). Since <p p is positively homogeneous of order 1, then 0* is an indicator 
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function: 



4>* p (r d \R d> - v; p) = sup (r d ' : i p ' + R d 'C - 4> p (v' ; V ; p 

(eP',C')er 

sup ( dev \T d '^ 



d*v(c:(e-e*)) ^ + Rd ,£_ {1 _ D) 



\dev(C:(£ 



2G g(v,p) //3 



2Gg(y,p)/l3 



-c 



(38) 



= I E (t* i^jvjp 

The indicator function 1^ is associated with the convex set 

E = |(V d ', such that / p (r d ' , R d '; v; p) < o} (see Figure 4b) with 



f p {r d \R d ' ] w-p)=dev (r d ') : 



,x dew (C : (e - e p )) 



:i-d) 



\dev (C : (e - £ p ) 



+ ,R d ' (39) 



2Gg(v,p)/P v y 2Gg(v,p)/P 

The function / p is the loading function for an endochronic model with plastic incompressibility 
and with isotropic damage. It is associated with the loading domain E. If the past-history 
parameter p is a scalar equal to e p , the plastic dissipated energy, then a work-hardening 
behavior is defined, in the sense that the loading function evolves with the plastic dissipated 
energy. A differe nt approach to define work-hardening plasticity models was proposed by 



Ristinmaal (119991 ). 



The generalized normality conditions imposed on <ft* leads to: 



.p _ ; df p (T d ,R d ;v;p) _ ■ dev{C:(e-eP)) 
S - A d T i' - A 2Gg{y,p)/p 



C = A 



0/ p (r d ,i? d ;v;p) 



A 



(40) 



\f p (r d ,R d ;v;p) =0 f p (r d , R d ; v; p ) < 0, A>0 



where the last three inequalities are the Kuhn- Tucker conditions. The plastic flow defined in 
(7) is retrieved. Note that the derivatives are taken with respect to the generic variables r d ' 
and R d \ but they are computed at the present state r d ' = r d and R d '=R d . In summary, the 
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usual notions of plastic multiplier and loading surface have been defined for an endochronic 
model with damage. This kind of thermo dynamic formulation for e ndochronic models is quite 



innovative and has been first presented in 



Erlicher and Point 



(120061 ). for the case of no damage. 



As was pointed out in that paper, an important property characterizing endochronic models 
is the fact that at the actual state, the loading function f p is always zero: for this reason, 
the consistency condition f p = is always fulfilled as an identity and cannot be used to 
compute the plastic multiplier A. This is also true in this case, where the actual state is 
(r d , R d ) = ((1 - D) C : (e - e p ) , 0) . As a result, the Kuhn- Tucker conditions reduce to 
A = ( > 0, where C is the flow of the internal variable associated with R d and, using the 
language of the endochronic theory, is also the flow of the intrinsic time measure; it can be 
freely defined, provided that it is non-negative and that rate-independence is guaranteed. As 
already observed, the standard choice is ( = \\dev 



Figure 5 illustrates an example of uniaxial behavior of an endochronic plasticity model with 
damage. The parameters of the elastic phase and of damage (Definition 1) are the same as 
those of Figure 3. In addition, g = 1, ( is given by ( 9) with n = 5 , (3 = 2834.9 MPa l ~ n 
and 7//? = -0.5 ; as a result, a u = (2G/ {(3 + 7)) 1/n = 2.25 * ^2/3 = 1.8371 MPa, where a u 
is the upper limit of \\dev (r d / (1 - £>)) || = \\dev (er/ (1 - £>)) || = ^2/3 <r n / (1 - D) when 
<? = !■ 



In the example of Figure 6, the damage is defined by the rule D = 1 — 1/(1 + c^e v \ with 
c v = 1500 m 3 MJ~ 1 (Definition 2). The parameter c v indicates the sensitivity of dam- 
age to the energy e p dissipated by plasticity. If c v is large, the damage increment at a 
given e p - value is larger than in the case of small c v . The Young modulus and the Pois- 
son ratio are the same as in the previous figures. The parameters defining the intrinsic 
time flow (9) are: n = 15, (3 = 16.1846 MPa 1 ~ n and 7//? = —0.8 ; as a result, a u = 
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(2G/ {(3 + 7)) 1/n = 2.25 * ^2/3 = 1.8371 MPa. Moreover, the hardening function is defined 
1 + \\dev (e (t')) || ma a; t / 6 r ot ]/ £ u) j where t is the present time and e u = 0.0002. Figure 
6d depicts the evolution of Y d = 1/2 (e — e p ) : C : (e — e p ), i.e. the actual value of Y d> . Ac- 
cording to (35), this quantity is also equal to l^' ax , which is the upper limit of Y d ' . The typical 
endochronic behavior with plastic strains increasing during unloading phases is highlighted in 
Figure 6b. As a result, owing to the damage rule depending on the dissipated plastic energy, 
also the damage slightly increases during the unloading phases: observe the damage evolution 
after t=0.3 and t=0.5, which are the instants where unloading phases begin. 



A brief discussion about stability and uniqueness 



It is well-known that standar d endochronic models violate the Drucker's postulate and the 



Ilyushin's postulate, see e.g. ((Sandler 



19781 ) . As a result, inelastic strains may continuously 
increase if a cyclic stress of constant and arbitrarily small amplitude is imposed around a given 
static stress (Figure 7a). Dually, a stress relaxation occurs when cycling straining is imposed 
(Figure 8a). The parameters used for the numerical simulations of Figures 7 and 8 are: E = 
35000 MPa, v = 0.18, g = 1, j/@ = -0.5, while (3 has a value such that (2G/ (f3 + 7)) 1/n = 
1.8371 MPa, for the given n values used in the figures. The strain accumulation entails a 
violation of a Lyapunov-type stability condition . For thi s reaso n, endochronic theory have been 



repeatedly criticized in the past years. However. 



Bazant 



(119781 . p. 705) showed that endochronic 



models do fulfil some weaker physically motivated stability conditions. Moreover, there are 
materials that are stable in the Drucker's sense and others that are not. Hence, for these 
materials, a proper model cannot fulfil the postulate of Drucker. All the aspects concerning this 



subje ct have been explored in detail in the previously cited references (ISandler 



1978 



Bazant 



19781 ) for endochronic models without damage. A detailed analysis for the case of models 
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with damage would deserve further studies, but this is beyond the purposes of this paper. 
Figures 7a and 8a simply show the influence of the parameter n on the strain accumulation 
and the stress relaxation for an endochronic model without damage. When n tends to infinity, 
a plastic behavior of Prandtl-Reuss type is retrieved, where neither strain accumulation nor 
stress relaxation occur. Figures 7b and 8b concern models with damage. 



Another important topic concerning p 



ness due to strain-softening; see e.g. (jJirasek and Bazant 



asticity and/or damage m odels is the loss of unique- 



2002). An exhaustive treatment 



of this subject for endochronic models with damage requires further analyses. However, 
for illustrative purposes, a simple analytical study of a uniaxial model is presented here- 
after. Let a, e and e p be the stress, the total strain and the plastic strain in the axial di- 
rection, respectively. Then, the uniaxial behavior can be represented by the following law: 
a = (1 — D) E (e — e p ) = (1 — D) Ee e , where E is the Young modulus. The incremental form 
reads 

da = (1 - D) E ids - de p ) - a-^- = (1 - D) Eds - (3a^ - a-^— (41) 

1 — D g 1 — D 



where the intrinsic time increment is 



d( = ^1 + -jjSign (ade)^j \a\ n 1 (1 — D) 1 n \de\ 



(42) 



and the damage increment writes 

dD = H{M f^j^ r) (43) 

with R(e e ) = E(e e ) 2 /2 and f D = {l-D) s R{e e ) - r < 0. Assume a > and de > 
(loading); the case a < 0, de < is analogous. Then, the condition to avoid strain-softening 
is 

da > (44) 
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The generic damage increment when fn = is given by 



2 1 — D 2 n — D) 2 

rfD = de e = -E- '- ids - de p ) 

s e e s a v ' 



Moreover, from (41) one has de p = /3ad(/ [E (1 — D) g] . Hence, the condition (44) assumes 
the following form 



1 — D) Eds — (3a—) (l — -) > (45) 



9 ) \ s 

The first factor is always positive provided that g > 1 . This can be proven using the definition 
of d( given in (42) with a > 0, de > and observing that the non-negativity of the first factor 
in (45) is equivalent to the condition a j (1 — D) < (E/ (f3 + 7))^™ (g) 1 ^ = a y (g) 1 ^, stating 
that the effective stress is always less or equal than the bounding axial stress a y , modified 
by the hardening function g. If g > 1, this inequality is always strictly fulfilled. Hence, 
strain-softening can be avoided if s > 2 . The same result can be obtained using the tensor 
expressions (7), (9), (27) and (33) and imposing that all the stress components are zero except 
(T11 : = °- This proof is omitted for brevity. 



case of elasticity with damage (INedjar 



?he same condition on s has been found for the 
200ll ) . Note that g < 1 induces strain-softening also 
when there is no damage. The analysis of the unloading case is not necessary, since at a given 
stress-strain state with a 7^ 0, the unloading stiffness is always greater than the loading one. 
A more complex analysis, not considered here, is needed for the multi-axial case, where the 
fourth-order tensor of tangential moduli for the endochronic model with damage should be 
computed. If strain-softening is avoided, the uniaxial behavior in what concerns the strain 
accumulation and the stress relaxation is analogous to that of standard endochronic models. 
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Conclusions 



An extended endochronic theory with a scalar damage variable was developed, based on the 
postulate of strain equivalence and by using pseudo-potentials depending on state variables 
and on parameters related to the past history of the material. The relevant loading surfaces, 
for damage and for plasticity, were defined. Two different damage pseudo-potentials were 
discussed and a formalization of the conditions on state variables affecting the definition of 
damage was provi ded, by an additional ind icator function in the Helmholtz free energy. In a 



companion paper (jErlicher and Bursil . 



20071 ) , a link between this extended endochronic theory 
and the Bouc-Wen type models with both strength and stiffness degradation is established. 
This will permit to prove the thermodynamic admissibility of these Bouc-Wen models and to 



highlight a constraint for the relevant stiffness degradation rules. 



Appendix: Notations 



The following symbols are used in this paper: 

C = fourth-order elasticity tensor 

D = internal variable associated with isotropic damage 

en — energy per unit volume dissipated through damage 

e p = energy per unit volume dissipated through plasticity 

fo = loading function for damage 

f p = loading function for plasticity 
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G = shear modulus 

g = hardening-softening function 

H(-) = Heaviside function 

I = fourth-order identity tensor 

In = indicator function of the set H 

K = bulk modulus 

q d = dissipative thermodynamic forces vector 

q rad = non-dissipative thermodynamic forces vector 

R d = dissipative part of the thermodynamic force introducing isotropic hardening(softening) 
R nd = non-dissipative part of the thermodynamic force introducing isotropic hardening(softening) 
v = state variables vector 

Y d = dissipative part of the thermodynamic force dual to the damage variable 
ynd _ non-dissipative p ar t of the thermodynamic force dual to the damage variable 
z = hysteretic part of the stress tensor 

f3 = coefficient defining the plastic flow of Endochronic models 
7 = coefficient defining the plastic flow of Endochronic models 
e = total small strain tensor 
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£ p = plastic small strain tensor 

£ = intrinsic time measure for Endochronic models. Moreover, it is the internal variable 
associated with isotropic hardening/ softening of Endochronic models 

d = intrinsic time scale for Endochronic models 

A = plastic multiplier 

\d = damage multiplier 

fi = hereditary kernel 

p = history-dependent parameters vector 

cr = Cauchy stress tensor 

cr d = dissipative part of the Cauchy stress tensor 
cr nd = non-dissipative part of the Cauchy stress tensor 

T d = dissipative part of the thermodynamic force dual to the plastic strain tensor 

T nd = non-dissipative part of the thermodynamic force dual to the plastic strain tensor 

$! = mechanical or intrinsic dissipation 

4> = pseudo-potential or dissipation potential for plasticity 

0* = dual pseudo-potential for plasticity 

4>d = pseudo-potential or dissipation potential for damage 

(f)* D = dual pseudo-potential for damage 
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^ = Helmholtz free energy volume density 
1 = second-order identity tensor 
()=McCauley brackets 
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List of Figures 



Classification of Endochronic models and most relevant references. The 
black square indicates that a thermodynamic formulation based on a suited 
pseudo-potential was found for the associated group of models. [34 

The increment de p of the energy dissipated by plastic strain (1) and the 
increment dep of the energy dissipated by damage (2). The increments of 
the elastic and plastic strain as well as of the stress are also schematically 
illustrated for a) the loading phase; and b) the unloading phase. [35 

Elastic with damage model. Uniaxial behavior, (a) The stress-strain evolution. 
The dotted line indicates the effective stress o"n/(l — D). (b) The damage 
evolution, with the damage increasing only during loading phases, (c) 
Time-evolution of Y d , i.e. the actual value of the thermodynamic force Y d ' , 
and of y^a X , i.e. the upper bound of Y d ' defined in (29). [36 

Endochronic model with damage in the uniaxial stress regime, (a) Several 
configurations of the effective domain D of the pseudo-potential 4> p . (b) 
Corresponding configurations of the loading domain E, associated with the 
loading function f p . The actual stress point always occurs at R d ' = R d = 
and, at this point, the flows e v and ( are defined by the normality condition. |37 

Uniaxial behavior of an endochronic plasticity model with damage (Definition 

1). (a) The stress-strain evolution. The dotted line indicates the effective 

stress o"n/(l — D). (b) Damage evolution, (c) See Figure 3c. [38 
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Uniaxial behavior of an endochronic plasticity model with damage (Definition 
2). (a) The stress-strain evolution. The dotted line indicates the effective 
stress 011/(1 — D). (b) The stress as function of the plastic strain. Observe 
that the plastic strain increases during unloading phases, (c) The damage 
evolution, with damage slightly increasing also during unloading phases, (d) 
Evolution of Y d vs. time. 

Strain accumulation, uniaxial behavior with stress varying between 1.05 and 
1.24 MPa. (a) Endochronic model without damage, with the intrinsic time 
(4). (b) Endochronic model with damage, with the intrinsic time (9) and the 
damage evolution given by (27) and (33), with tq = 0.000012 and s = 6. 

Stress relaxation, uniaxial behavior with strain varying between 0.000075 and 
0.00008. (a) Endochronic model without damage, (b) Endochronic model with 
damage, with the same damage parameters as in Figure 7b. 
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Fig. 1. Classification of Endochronic models and most relevant references. The black square indicates 
that a thermodynamic formulation based on a suited pseudo-potential was found for the associated 
group of models. 
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Fig. 2. The increment de p of the energy dissipated by plastic strain (1) and the increment deo of 
the energy dissipated by damage (2). The increments of the elastic and plastic strain as well as of 
the stress are also schematically illustrated for a) the loading phase; and b) the unloading phase. 
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Fig. 3. Elastic with damage model. Uniaxial behavior, (a) The stress-strain evolution. The dotted 
line indicates the effective stress o"n/(l — D). (b) The damage evolution, with the damage increasing 
only during loading phases, (c) Time-evolution of Y d , i.e. the actual value of the thermodynamic 
force Y d ' , and of Y^ ax , i.e. the upper bound of Y d ' defined in (29). 
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Fig. 4. Endochronic model with damage in the uniaxial stress regime, (a) Several configurations of 
the effective domain D of the pseudo-potential <p p . (b) Corresponding configurations of the load- 
ing domain E, associated with the loading function f p . The actual stress point always occurs at 
R d ' = R d = and, at this point, the flows eP and ( are defined by the normality condition. 
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Fig. 5. Uniaxial behavior of an endochronic plasticity model with damage (Definition 1). (a) The 
stress-strain evolution. The dotted line indicates the effective stress <r 11 /(l — D). (b) Damage evolu- 
tion, (c) See Figure 3c. 
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Fig. 6. Uniaxial behavior of an endochronic plasticity model with damage (Definition 2). (a) The 
stress-strain evolution. The dotted line indicates the effective stress — D). (b) The stress as 

function of the plastic strain. Observe that the plastic strain increases during unloading phases, (c) 
The damage evolution, with damage slightly increasing also during unloading phases, (d) Evolution 
of Y d vs. time. 
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Fig. 7. Strain accumulation, uniaxial behavior with stress varying between 1.05 and 1.24 MPa. (a) 
Endochronic model without damage, with the intrinsic time (4). (b) Endochronic model with damage, 
with the intrinsic time (9) and the damage evolution given by (27) and (33), with tq = 0.000012 and 
s = 6. 



40 



(a) (b) 




Fig. 8. Stress relaxation, uniaxial behavior with strain varying between 0.000075 and 0.00008. (a) 
Endochronic model without damage, (b) Endochronic model with damage, with the same damage 
parameters as in Figure 7b. 
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